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Abstract. Beta Beam, the concept of generating a pure and intense (anti) neutrino beam by letting accelerated radioactive ions 
beta decay in a storage ring, called Decay Ring (DR), is the base of one of the proposed next generation neutrino oscillation 
facilities, necessary for a complete study of the neutrino oscillation parameter space. Sensitivities of the unknown neutrino 
oscillation parameters depend on the Decay Ring's ion intensity and of it's duty factor (the filled ratio of the ring). Therefore 
efficient ion production, stripping, bunching, acceleration and storing are crucial sub-projects under study and development 
within the Beta Beam collaboration. Specifically the feasibility of these tasks as parts of a Beta Beam implementation at CERN 
will be discussed in this report. The positive impact of the large #13 indications from T2K on the Beta Beam performance will 
also be discussed. 
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INTRODUCTION 

The discovery of neutrino oscillations has confirmed 
that neutrinos are massive and that their flavor (v e , v^, 
V T ) and mass eigenstates (Vi, V2, V3) are mixed. The 
first oscillation evidence of atmospheric muon neutrinos 
into tau neutrinos was made by Super-Kamiokande in 
1998 [1]. Neutrino physics is now in an era of preci- 
sion measurements of the parameters that govern these 
oscillations: two Am?- = m 2 . — m\. parameters (lAm^J, 

lAm^J), three mixing angles (#23, #12, #13) and a CP 
violating phase (Sep). The oscillation phenomena natu- 
rally fall into two domains: atmospheric, — )> V T , and 
solar neutrino oscillations, v e — )> V)^ T . In the solar sec- 
tor the latest la best fits are 2 Ara^ = (7.59+o; 2 g) x 
10- 5 eV 2 and G l2 = (34.0+^)° and in the atmo- 
spheric sector |Am^| = (2.50+°;^) x 1(T 3 eV 2 and 

#23 = (46.II402) [2]. According to recent indications 
from the T2K experiment #i 3 > 5° with 2.5a [3]. A 
global fit with the latest MINOS results suggests that 
1.8° < 0i3 < 10.8° with 3a [2]. If it will be proven that 
013 > it opens the possibility of discover 8cp ^ 0° and 
7^ 180° which would mean the existence of CP viola- 
tion in the leptonic sector. Another important unknown 
value is the sign of Am 2 : . A complete study of all neu- 
trino parameters requires a better characterized neutrino 
beam with higher flux than ever available before. One 
of the three present options for the next generation neu- 
trino oscillation facility [4] is the Beta Beam concept [5] 
wherein it is proposed to store high energy (7 = 100) ra- 
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2 For these values normal hierarchy, sig^Ara^) > 0, is assumed. 
See [2] for a complete set of values. 



dioactive ions in a racetrack shaped storage ring with a 
straight section pointing to a neutrino detector. The ions 

(a) 




FIGURE 1. Beta Beam baseline implementation at CERN 
for low Q ions ( 6 He and 18 Ne) (a) and high Q ions ( 8 Li and 
8 B) (b). 

have to have optimized decay time, to maximize the neu- 



trino flux, and Q-value 3 , for the chosen baseline, and to 
produce both a neutrino and anti-neutrino beam a pair of 
j3~ and j3 + decaying ions have to be used. Ions that beta- 
decay in the straight section emit electron (anti) neutri- 
nos in a pure (almost only v e or v e ) beam with an open- 
ing angle of 1 jy. The aimed annual (anti) neutrino flux 
of (2.9el8) l.lel8 from (j3 _ ) j3 + decaying ( 6 He) 18 Ne 
ions gives Sep sensitivities shown in fig. 2. The distance 
to the first v e —> -oscillation maxima for this ion-pair 
(with Q-values around 3 MeV) and thereby to the lo- 
cation of the detector is about 130 km. Corresponding 
sensitivity studies for an ion-pair with higher Q-values 
(around 13 MeV), 8 Li and 8 B, have been made and show 
that about a factor 5 times higher neutrino flux is nec- 
essary for similar sensitivities [6]. The baseline for this 
case is around 700 km. For both the low and high Q- 
value ion pairs the baselines studied in this report start 
at CERN and use PS and SPS as part of the accelerating 
scheme (see fig. 1). The sensitivities depend also on the 
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FIGURE 2. The 13 (left) and 8 CP (right) sensitivities for 
different suppression factors (SF) [6]. A SF of 0.01% (green) 
gives almost the same result as the case of no atmospheric back- 
ground (orange). The 0.1% SF (blue) gives similar sensitivity 
and with 1 % SF (red) the sensitivity has decreased slightly. 

suppression factor (SF) of the detector, which is the same 
as the duty factor for the storage ring. This storage ring is 
called the "Decay Ring" (DR) and has the same circum- 
ference as SPS, C = 6911.6 m. The required SF and its 
impact will be discussed in this report after ion produc- 
tion, preparation, bunching and acceleration is presented. 



different production schemes depending on the ion. For 
the 6 He ion production the ISOLDE [7] method will be 
used. High energy protons will bombard a heavy spal- 
lation target (either tungsten or lead) so that spallation 
neutrons are forced out to and react with the surrounding 
ISOL target (normally BeO) so that 6 He ions are pro- 
duced (see figure 3). This is a standard method and the 
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FIGURE 3. Spallation target (W or Pb) and ISOL target 
(BeO) for 6 He production. 

production rate can be estimated by scaling and would 
depend only on the proton source. Assuming a 200 kW 
Super Proton Linac (SPL) that gives 2 GeV protons or 
a 20 kW Rapid Cycling Synchrotron (RCS) 4 , also giv- 
ing 2 GeV protons, the production rate would be 5 x 10 13 
6 He/s. 

The 18 Ne production would need 160 MeV protons 
from a linac, possibly from an upgraded Linac4, to hit 
a molten salt (NaF) loop. Two reactions would then 
contribute to the 18 Ne production; l9 F(p,2n) ls Ne and 
23 Na(p,X) m Ne (see figure 4). According to studies 



2l/s 



7.5x25x1 5cm 




Protons 



6mA, 160MeV 



18Ne, 10 13 /s 



Temperature 
regulated 
transfer Line 
to ion source 



Irradiation 
cell 

800 kW 







[Dump 




200 kW 









AT=100°C 







Heat 




Diffusion 




Exchanger 


Chamber 




1 MW 





40x1 5x1 5cm 

e>90% 

Psalt = 10mBar 



NaF/ZrF4 

5 C 



ION PRODUCTION 

The CERN implementation of the baseline for the ion 
pair with low Q-value, 6 He and 18 Ne, will consist of two 



FIGURE 4. The molten salt loop to produce i8 Ne with two 
reactions; l9 F(p,2n) ls Ne and 23 Na(p,X) l ^Ne. If an upgrade 
of CERN's Linac4 gives 160 MeV protons this target station 
could give 10 13 18 Ne/s [9]. 



3 The ion's reaction energy. It defines the neutrino energy. 



4 This RCS are under consideration for the LHC upgrade and Hie- 
ISOLDe [8] would also benefit from it. 



1 MW target station (only twice as powerful as the cur- 
rent operating target station for CNGS) would provide 
about 10 13 18 Ne/s. This is however a novel idea and a 
proof of principle experiment is scheduled for Novem- 
ber at CERN [9] . 

Another novel concept is the production of 8 Li and 
8 B ions with the so called "Production Ring" (PR) [10]. 
For the 8 Li production 25 MeV 7 Li ions enters the PR 
that is a small storage ring with a wedge shaped gas-jet 
deuterium target (see figure 5). Multiple passes through 
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FIGURE 5. Production ring for 8 Li and 8 B production [11]. 

this target cause ion-cooling, stripping and production 
through the reaction ^H+3Li^ 8 Li+/?. The produced 8 Li 
ions would then exit the target with an angular distribu- 
tion and be collected by 2 ]Um thick catcher foils and then 
diffused away with a cold finger (see figure 6). The effi- 




FIGURE 6. Collection device with thin, 2 /im, catcher 
foils [12]. 

ciency of this collection device are to be determined in 
ongoing studies [12] . The same concept would be used 
for 8 B production but with the reaction ^He+^Li— ^ 8 B+?r. 
The cross-section and angular distribution of this reac- 
tion is however unknown. Therefore measurements have 
been performed and analysis is ongoing [13]. Produc- 



tion Ring studies have shown that there are difficulties to 
reach the required target density [14] therefore new ideas 
that suggests direct kinematics through a 7 Li film target 
will be tested [15]. Further details about ion production 
for the Beta Beams will be shown in [16]. 

BEAM PREPARATION 

The gas of the radioactive atoms from the ion production 
(see previous section) diffuse into a volume with a mag- 
netic field, B. This field makes free cold electrons move 
in a circular motion with the cyclotron angular frequency 
co ce = eB/m where e is the elementary charge and m is 
the mass of the electron. An electro magnetic wave with 
the same high frequency, CQhf = &>ce> is applied and al- 
low Electron Cyclotron Resonance (ECR) heating. Due 
to their subsequent energy, when they collide the atoms 
of the gas, they are able to ionize them. This ECR ion 
source has best efficiency producing singly charged ions 
so a stripper giving fully stripped ions is assumed after 
the ECR. R&D is ongoing for a 60 GHz ECR source 
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FIGURE 7. Magnetic field and temperature map of the 
60 GHz ECR source [17]. 

(see figure 7), which would be the most powerful ECR 
source to date [17]. It operates in pulse mode so that 
50 /is pulses of monocharged ions will be extracted from 
the ECR. The estimated efficiencies for the Beta Beam 
ions in this case are 30% for 6 He, 20% for 18 Ne and 3 to 
10 times less for 8 Li and 8 B. 



BUNCHING AND ACCELERATION 

Fully stripped ions enter the Beta Beam acceleration 
complex in 50 /is dc bunches and with about 8 keV/u. 
Prebunching and preacceleration is done with a Radio 
Frequency Quadrupole (RFQ) and a Linac. A moderate 
acceleration gradient of 3-6 MV/m is assumed so nor- 
mal RF cavities can be used. The Linac will accelerate 



the ions up to 100 MeV/u before they are transferee! to a 
Rapid Cycling Synchrotron (RCS). The first task for the 
RCS is to create bunches with high ion intensity. This 
is done by multiturn injection. With the current design 
of the RCS for the Beta Beams [18] the revolution time 
is 1.96 /is. The ECR pulses are 50 /is so the injection 
takes 26 turns. Thanks to phase- space rotation the 26 in- 
jections fit into one single RCS bunch with 50 % effi- 
ciency. The second task of the RCS is acceleration up to 
about 15 Tm before the high intensity bunch is transfered 
to CERN's Proton Synchrotron (PS). The PS collects 20 
RCS bunches one by one. It will be 1.9 seconds between 
the first and last injection. There will therefore be inten- 
sity differences in the bunches when all 20 bunches are 
injected to the PS due to radioactive decay (see figure 8). 
After injection the batch of 20 bunches are accelerated up 
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FIGURE 8. Single bunch intensities are indicated in green 
for the 1st, 5th, 10th and 20th bunch and the total number of 
helium ions are indicted in red. The schematic time structure 
of the RCS, PS and SPS magnet cycles are also shown in 
black [19]. 

to maximum PS capacity, 86.7 Tm, and then transfered 
to CERN's Super Proton Synchrotron (SPS). For proton 
operation the PS bunches are normally split before enter- 
ing SPS so that they can be handled by the 200 MHz (h 
= 4620) SPS RF system. However, for the Beta Beams 
we assume that an additional slower, 40 MHz (h = 924), 
RF system has been implemented into the SPS. This 
is to mitigate space charge effects. Beta Beams have a 
high ion intensity and specially a high charge intensity. 
For example for nominal Beta Beam neutrino fluxes we 
need 4.9xl0 12 charges per bunch in the SPS and the 
comparable most intense beam accelerated in the SPS 
is 3.5 x 10 11 charges per bunch. Space charge effects are 
worse for low energies (it goes as y -2 [20]). Therefore 
large bunches (same size as given by PS) are taken care 
of by the slow 40 MHz and 1 MV RF system at SPS 
injection where 7=9. After initial acceleration space 
charge effect is less crucial and the bunches are short- 
ened so that the RF can swap to the normal 200 MHz and 
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FIGURE 9. A simplified diagram showing how one acceler- 
ation part of the SPS RF cycle could look like. 



8 MV RF system. With this RF system all ions can then 
be accelerated across transition and then up to a y = 100, 
for all ions. 



ACCELERATED INTENSITIES 

The ion intensities through the whole Beta Beam com- 
plex are calculated [21] starting with the source rate and 
taking into account the radioactive decay and the effi- 
ciencies of all components (discussed above). Results 
are shown in table 1 . They are shown for the low-Q ion 
pair only since source-rate and efficiencies are still too 
uncertain for the high-Q ion pair. Intensities are shown 
per bunch at the end of each components cycle, i.e. the 
last row shows the intensities in SPS when the ions have 
reached y = 100. The injection into the DR and the DR 
intensities will be discussed in the next section. 
TABLE 1. Ion bunch intensities after each stage in the Beta 
Beam complex before the Decay Ring. Corresponding values 
for 8 B and 8 Li are not shown here due to large uncertainty in 
source rate and ECR efficiency. 
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DECAY RING 

The Decay Ring (DR) is a racetrack shaped storage ring 
aiming the neutrino beam from the decaying ions to a 
neutrino detector in a big cave somewhere in Europe. 
For the low Q ions the DR will have a 0.6° declination 
angle pointing towards the Frejus tunnel (French/Italian 
border) and for the low Q ions the DR will decline with 
a 2.9° or 3.3° angle pointing towards Canfranc (Spain) 
or Gran Sasso (Italy) respectively. In the current de- 
sign [22] the DR has the same circumference as SPS, 
6911.6 m, and the efficient straight section is 39%. This 
gives bending radii as small as 121 m so that supercon- 



ducting dipole magnets, giving a 4 to 8 T magnetic field, 
are necessary. One important aspect of the DR is that the 
ion bunches can only fill a small part of it's total length. 
Before the hints of large #13 (see next section) it was as- 
sumed a duty factor as small as 0.58% is necessary for 
physics reach. The 20 bunches then had to be limited to 
a bunch length of 691 1.6m x 0.5 8%/20 = 2 m each. There 
is also a strong requirement of ion intensities in the DR. 
With 20 bunches, the number 18 Ne ( 6 He) per bunch have 
to be 3.4 x 10 12 (4.5 x 10 12 ) (and about 5 times more of 
8 B ( 8 Li) ions) to reach the nominal (anti) neutrino fluxes. 



Decay Ring Injection 

A dedicated injection scheme of the SPS bunches into 
the DR has been developed [23, 24] in order to increase 
the ion bunch intensities and limit the bunch lengths. 
This scheme uses two RF systems and includes catch- 
ing the SPS bunch injected off-momentum, merging the 
fresh SPS bunch into the existing DR bunch (to increase 
the bunch intensity) and collimation at dp/p = 2.5%o 
(to limit the bunch length to 2 m). Numbers of ions 
per bunch injected from SPS are shown in the last row 
in table 1. During the first injections only ions that are 
not captured are collimated away (capture efficiency is 
about 90% for all ions [24]). Due to RF-gymnastics dur- 
ing merging the phase- space of the DR bunch increases 
and after about 15 merges many ions are collimated af- 
ter each injection so that the bunch length do not exceeds 
2 m. Due to this collimation and to the constant decay 
of the radioactive ions the number of ions per bunch 
saturates after about 20 injections. For the 6 He case 
where 3.8 x 10 11 ions are injected per bunch the accumu- 
lation of the DR bunch saturates at about Nu at = 3 x 10 12 
6 He/bunch i.e. below the required 4.5 xlO 12 6 He/bunch 
(see figure 10(a)). From the number saturated ions per 
bunch we get the annual neutrino flux from 
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where T e ff = 10 7 seconds 5 and all other parameters are 
given in table 4 and 3. The results from these studies 
and for SF = 0.58% are summarized in first and second 
column in table 4 for 18 Ne and 6 He respectively. 



Decay Ring Storage 

High intensity bunches can have non-negligible 
amount of charges which could cause the particles to 
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FIGURE 10. (a)Accumulation of 6 He ions in a DR bunch. 
Saturation occurs at 3x 10 12 6 He/bunch, less than the required 
4.5 xlO 12 6 He/bunch. (b) Bunch intensity limits according to 
HEADTAIL [25], MOSES [26] and the Coasting Beam equa- 
tion [27] for 18 Ne in the DR. 



interact with each other and with the vacuum chamber. 
These, so called "Collective Effects", could cause many 
different reasons for beam instability. For example "sin- 
gle bunch resonance broad band impedance" appears 
when electro magnetic wake fields from the head of 
the bunch are captured by cavities in the vacuum pipe 
and interact with the t ail of the bunch. If the quality 
factor is Q = Ry/C/L and the resonance frequency is 
(O r = 1 / ^/LC the resonance impedance can be modelled 
as an RLC circuit [28] in the transverse plane as 



where the transverse shunt impedance, R±, is a value 
indicating the total divergence from a perfectly smooth 
vacuum pipe around the whole ring. The value for RHIC 
is about R R ^ IC « 2 M^/m [29]. The DR will be a mod- 
ern machine so we will assume it will have a smooth 
vacuum pipe design and therefore twice as good shunt 
impedance; R^ R w 1 M^/m [30]. However, simula- 
tions, modeled on eq. (2), with a tracking program called 
HEADTAIL [25], theoretical calculations, with a pro- 
gram called MOSES [26], and estimations with a sim- 
plified equation called Coasting Beam equation (CB 
eq.) [27] all claim that although « 1 M£l/m reso- 
nance impedance would cause instable beams below the 
required 3.4 x 10 12 18 Ne/bunch (see figure 10(b)). For 
example, according to HEADTAIL the 18 Ne ion beam 
would be unstable if there were more than 10 12 ions per 
bunch, i.e. more than a factor 3 below the required. All 
parameters used for these simulations and calculations 
are listed in table 2 and 3. 



TABLE 2. Input parameters (some from [31]). Parame- 
ters below the line are all calculated. 



Parameters Description 


DR 18 Ne DR 6 He 


Z Charge Num. 
A Mass Num. 
h/2 M Half Life 
Vrf [MV] Voltage 
E rest [MeV] Rest Energy 
Lb [m] Bunch Length 
8 max [10 -3 ] Mom. Spread 
e^(lcr) [7rmxrad] 
£n (la) [7rmxrad] 


10 2 
18 6 
1.67 0.81 
26.75 
16767.10 5605.54 
1.97 1.97 
2.5 2.5 
1.48e-05 1.48e-05 
7.90e-06 7.90e-06 


E tot [GeV] = yxE r£?s , 

% [ns] = f c 

as[e- 3 ] = 5 -T 

e] a [eVs] = jP 2 E tot T b 8 max 

&[io-3]=y^« 

r [am] = r p Z 2 /A 


1676.7 560.6 
6.57 6.57 

1.25 1.25 
43.20 13.95 
8.1 3.7 
8.53 1.02 



IMPACT OF LARGE 13 

In the previous section we showed two reasons why the 
DR could not use the ions provided by the rest of the 
Beta Beam complex efficiently enough to reach required 
neutrino flux; Not enough ions could be accumulated per 
bunch with the merging technique and the bunch inten- 
sity limit due to collective effects were also below the re- 
quired. The underlying reason for both of these effects is 
the limited suppression factor; 0.58%. New results from 
the T2K experiment show however indications [3] that 
would relax the duty factor requirements in the DR and 



TABLE 3. Input parameters (some from the previous 
Beta Beam Decay Ring design report [31] and some 
from updated design [32]) above the first line. Assumed 
transversal impedance parameters between the lines. Cal- 
culated parameters below the last line. These parameters 
are the same for the different isotopes. 
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*eff 
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Ytr 
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18.6 


1 




Relativistic Gamma 


100.0 


Qx 
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21.23 


Qy 




Vertical Tune 


17.16 


{P)x [m] 




rXv. A-jJlfUrl r UrlC 


1 94 70 


(P)y [m] 




Av. y-f5tron Func. 


160.40 


{D) x [m] 




Av. Dispersion 


-0.60 






x,y Chromaticity 


0.0 


b x [cm] 




x Pipe Size 


16.0 


by [cm] 




y Pipe Size 


16.0 


Q± 




Quality Factor 


1.0 


fr [GHz] 




Resonance Frequency 


1.0 






Shunt Impedance 


1.0 


P = V^ 


y-2 


Relativistic Beta 


1.00 


*7 = Ytr 2 ~ 


y-2 


Phase Slip Factor 


2.8e-3 


T rev [ps]= 


C 
Jc 


Revolution Time 


23.06 


frevim = 


YZ, 


Revolution Frequency 


0.27e6 


R [m] =C/2k 


Machine Radius 


1100 


© c [GHz]= 


. Pc 

'Ky 


Cut-Off Angular Freq. 


1.87 



open new opportunities for the Beta Beams. Figure 2 
shows that in the #13 region allowed by T2K's 6 v e can- 
didate events (green zone indicated in fig. 2) the sup- 
pression factor is much less crucial. A suppression factor 
of 2% (black line) would give more or less as good Sep 
sensitivity as SF = 0.5%. The impact of filling four times 
more of the DR, 2.3% instead of 0.58%, has therefore 
been studied. This opportunity could be used in differ- 
ent ways but so far only one has been studied; to fill the 
DR with the same bunches, 2 m long, but with 4 times as 
many bunches. This simple solution can be achieved by 
letting SPS wait for 4 PS batches before the SPS batch of 
80 bunches is injected into the DR. More ions will then 
have decayed and average SPS bunch intensity will be 
lower; 0.9xl0 n 18 Ne/bunch and 2.5xl0 n 6 He/bunch. 
The average DR bunch intensity will therefore saturate 
at a lower value (compare figure 1 1(a) with figure 10(a)). 
The required bunch intensity has however decreased with 
a factor 4 (since we have 4 times more bunches) so that 
the accumulation succeeds to provide with enough neu- 
trino flux. The result for these studies and for SF = 2.3% 
are summarized in third (for 18 Ne) and fourth (for 6 He) 
column in table 4 where also comparison with the nomi- 
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TABLE 4. Capturing and merging simulation results for the 
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low-Q isotopes after injected bunch intensities on the 2 nd row 
are assumed. Neutrino fluxes are calculated using eq. (1). 
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periment will take place at CERN in November 2011. 
There are well established beam preparation studies, for 
example a state of the art studies for a 60 GHz ECR is on- 
going. The new indication of large 613 has opened many 
new possibilities for the Beta Beams. The result of the 
impact- studies of one of them were presented here. We 
saw that the DR now can have injected and store enough 
ions to reach required neutrino flux for the low-Q ion pair 
baseline. 
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FIGURE 11. In case of large i3 SF can be relaxed to 2.3% 
and then 6 He saturation (a) occurs at 1.4 xlO 12 6 He/bunch, 
more than the required 1.1 xlO 12 6 He/bunch. Bunch intensity 
limits due to collective effects (b) also relax above required. 



nal neutrino fluxes are made. 

Thanks to the decrease of required bunch intensity 
the 80 ion bunches can also store enough ions before 
resonance impedance cause beam instability (see fig- 
ure 11(b)). However other collective effects reasons will 
need to be studied to confirm this. 



CONCLUSIONS 

A detailed study program for Beta Beam implementa- 
tion at CERN has been presented in this report. There 
has been encouraging ion production progress towards 
required source rate, specifically the 18 Ne production ex- 
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